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ABSTRACT 

We report on our first set of spectroscopic Hubble Space Telescope observations of 
the z ^ 11 candidate galaxy strongly lensed by the MACS J0647.7+7015 galaxy cluster. 
The three lensed images are faint and we show that these early slitless grism observations 
are of sufficient depth to investigate whether this high-redshift candidate, identified by 
its strong photometric break at 1.5/im, could possibly be an emission line galaxy at 
a much lower redshift. While such an interloper would imply the existence of a rather 
peculiar object, we show here that such strong emission lines would clearly have been 
detected. Comparing realistic, two-dimensional simulations to these new observations 
we would expect the necessary emission lines to be detected at > 5a while we see no 
evidence for such lines in the dispersed data of any of the three lensed images. We 
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therefore exclude that this object could be a low redshift emission line interloper, which 
significantly increases the likelihood of this candidate being a bona fide z ^ 11 galaxy. 


1. Introduction 


The search for the first galaxies has been at the forefront of extragalactic astronomy ever since 
the initial discovery that QSOs are so distant that their UV spectral features were redshifted into 


the optical pass-bands (Schmidt 1963). Owing to their high space density, star forming galaxies 


comprise a large fraction of the high-z galaxy population available to us at the highest redshifts. 
Gravitational lensing is an efficient method that allows the detection of high redshift star forming 
galaxy candidates that would otherwise have luminosities that are too faint to be detected using 
current space- or ground-based observatories ( Wyithe et al.|[MTT ). 


A strong z ^ 11 candidate was identified by Coe et al. (2013) in observations of the 2 ; = 


0.591 galaxy cluster MACSJ0647.7-h7015. These observations were made as part of the Cluster 
Lensing And Supernova survey with Hubble (CLASH, [Postman et al. 2012) program. Analysis of 
three separate lensed images (JDl, JD2 and JD3) of the faint red object in near infrared pass- 
bands led Coe et al. ( 2013D to conclude that this target was a strong candidate for a 2 ; 11 

galaxy. Three images of the same high redshift candidate were identified with magnification factors 
estimated to be 8,7, and 2 for the JDl, JD2, and JD3 images, respectively. The three images 
are magnified to F160W AB magnitudes of ~ 25.9,26.1, and 27.3, respectively. All three images 
show no detected flux in observations taken in 15 different HST filters over a wavelength range of 
0.2-1.4/im. However, all three were detected in the F140W (> 6a at 1.4/im) and F160W (> 12a at 


1.6/im) filters. While Coe et al. (2013) ruled out that this source could be a low-redshift interloper 
by considering a variety of possible low-redshift objects and the strength of the lensing model, 
their conclusion relied on the detection of this object in only two infrared broad band filters, while 
this source remained undetected at the longer wavelengths probed by relatively low signal-to-noise 
Spitzer/IRAC observations at 3.6/im and 4.5/im. The authors pointed out that the discovery of 
such an object was consistent with expectations extrapolated from 2 ; ~ 8 ( [Bradley et al.|[^012| ) but 
is in conflict with the significantly lower than expected number of 2 ; > 9 detected by [Bouwens et 


al. (2014) and Oesch et al. (2013). Results from Bradley et al. (2012); Oesch et al. (2013) were 


indicative of a dramatic buildup in the number of galaxies and of the cosmic star formation rate 
density over a very short period of time (< 200 Myr) between z ^ W and 8. This was later 


repeatedly confirmed by subsequent observations as well as re-analyses of earlier data (e.g. Oesch 


et al.||2014a|b; 

Ishigaki et al.||2014; 

Zheng et al.||2014; Bouwens et al.||2014). There is currently only 

one other multiply imaged galaxy 
to date is the z ^ 12 candidate U 

at 2 ; 10 (Zitrin et al. 

2014).The only other 2 ; > 11 candidate 

DF12-3954-6284 ( 

Ellis et al. 

2013 

). This source has only been 


detected in a single band and it can therefore be more easily explained away as being a lower- 


redshift interloper with strong nebular emission (Pirzkal et al. 2013). However, UDF12-3954-6284 


remains a possible high-redshift candidate as no emission lines have been conclusively detected in 
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subsequent observations and analyses (Bouwens et al. 2013; Brammer et al. 2013; Capak et al. 


2013). 


Confirming the nature of high-redshift candidates requires spectroscopic observations. If some 
of the very few z > 9 candidates already identified were to be shown to be intermediate redshift 
sources, then the observed deficit of galaxies at z 10 would be even greater, requiring an even 
more rapid buildup of early galaxies. Unfortunately, spectroscopic confirmation of high-redshift 
candidates has lagged behind due to the faintness of the sources, the apparent lack of Ly-a emis¬ 
sion in high-redshift galaxies, and the bright infrared sky background in ground based observations 


troscopic confirmation with current telescopes and limiting the spectroscopic studies possible with 
future telescopes. Fortunately, cluster lensing magnifies images of the MACS0647-JD candidate 
such that space-based low-resolution spectroscopy is feasible. 


(Shibuya et al. ~| 2012[[Ono et al.]|2012 |. The unlensed Hubble Ultra Deep Field (HUDF [Beckwith] 


et al. |200^ z > 9 candidates are significantly fainter (29.3 < F160W AB < 29.7), prohibiting spec- 


While the possibility that this object is a low redshift emission line galaxy is in contradiction 
with current lens models, spectroscopically ruling out that this object is not a low redshift star 
forming galaxy is important because our treatment of the nebular emission might be simplistic and 
incomplete. There is empirical evidence for the existence of faint moderate redshift galaxies with 
very strong emission lines. Straughn et al. (2011) identified several extreme emission line objects 
in Early Release Science observations using the G102 and G141 WFC3 grisms. Their Figure 3, for 
example, shows the observed spectrum of a F098M 26.87 AB magnitude source (ID397) at z=1.76 
and which has a very strong [OIIIJ/H/3 doublet (9.69 ± 1.3 10“^^erg/s/cm^, observed 5000A). 

This emission line is sufficient to produce a photometric break of nearly two AB magnitudes and 
is nearly as extreme as what would be required to reproduce the flux measured in MACS0647-JD 
5 X 10“^^erg/s/cm^). See also Frye et al. (2012) for other examples of extreme emission line 

10“^^erg s“^ cm“^. 


( 

galaxies with large single emission line fluxes of up to 1.38 


Redshift confirmation is therefore essential to definitively rule out a low redshift emission line 
interloper (or perhaps discover the true nature of a very interesting galaxy). 


Observations 


This paper is based on a set of new Hubble Space Telescope near-infrared stilless spectroscopic 
observations of the MACSJ0647.7+7015 cluster. As part of a follow-up project, we will be observing 
this cluster for 12 orbits using the G141 grism on the Wide Field Camera 3 (WFC3). These 
observations are split in three epochs, revisiting this field at different position angle on the sky. 
The position angles were chosen to minimize the amount of spectral contamination in the spectra 
of our target and to provide stronger, independent detections of any spectral break or emission 
lines in the spectra of this source. We have obtained our first set of observations, which we 


present here. The data were processed using the aXe data reduction package (Pirzkal et al. 2001 
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Kiimmel et al.|[2009 ), using a custom sky subtraction step which includes different components for 


Zodiacal light, Earth limb, and stray light. More details on this sky subtraction can be found in 


Brammer et al. (2014). We used archival data to create mosaics of the field in the F105, F125W, 


F140W and F160W filters and to model the SED of every object in the field. This allowed us 
to model the spectral contamination caused by nearby objects using aXe’s FluxCube quantitative 
estimates. This approach also allowed us to perform both regular and optimally weighted spectral 
extractions. Estimating spectral contamination is crucial when analyzing slitless spectroscopic 
because overlapping spectra can cause features that look very similar to spectral breaks and zeroth 
dispersion orders of nearby sources can also look similar to strong emission lines. A complete 
understanding of the field is therefore needed to rule out spectral contamination as the cause of 
any feature we detect in a spectrum. 


3. High redshift or Low Redshift Emission Line Interloper? 


While Coe et al. (2013) suggest the most plausible explanation for the detection of MACS0647- 


JD in the F140W and F160W filters is that the object is a lensed galaxy at z=ll, it is not the only 
possible scenario. Amongst the possible scenarios considered by the authors, one consists of the 
coincidental presence of two emission lines, one in each filter, to mimic a photometric break. 


To further explore this, we started by re-fitting the data from Table 2 of Coe et al. (2013) using 


the Monte Carlo Markov Chain approach (MCMC) described extensively in Pirzkal et al. (2012) 
and used to analyze other high-redshift candidates ( [Pirzkal et al. 2013). Our approach uses Bruzual 


& Chariot (2003) stellar population models, a Salpeter initial mass function, and SSP models. The 


free parameters are the redshift, stellar mass, stellar age, extinction level (Calzetti et al. 2000), 
metallically and the ionizing photon escape fraction. This new analysis confirmed that the Coe 


et al. (2013) photometry were better fit, especially at bluer wavelengths, by a template consistent 
with a high-redshift Lyman break object dit z ^ 11. Figure [^demonstrates that a IO^^Mq galaxy 
at a redshift of z ^ 11 provides a very good fit. Although, as we show, such a model would require 
a significant amount of flux, about the same level that was detected in the F160W filter, in the 
IRAC bands at 3.6 and 4.5/im, the IRAC data are not very sensitive and are consistent with a 
high-redshift Lyman-Break galaxy model. 

Since the MCMC SED fitting approach allows us to fully explore our input parameter space, we 
can examine what low redshift models {z < 2.5) best fit the photometric observations. However, this 
approach is ultimately limited by our handling of the nebular emission process, which we assume 
here to be the source of the high Equivalent Width (EW) lines in the spectra of our simulated 
star forming galaxies. In Figure we show two distinct cases, one at z ^ 1.4 and one at z ^ 2.2, 
where a young (few 10^ year old) star forming galaxy with strong nebular activity can reproduce 
the broad band observations. In our treatment of nebular emission, we have assumed the same 
ionizing photon escape fraction as for the nebular emission. The detailed description of how nebular 


emission is simulated can be found in Pirzkal et al. (2012). Ultimately, it is the ratio of the nebular 
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line emission to the continuum emission which limits how much of the broad-band photometric 
break can be reproduced by the nebular emission lines of our models . We also must resort to 
an amount of extinction large enough (Ay ~ 3-4) to quench the amount of blue light, including 
the nebular continuum light, so that such an object can remain undetected in broad band filters 
blue-ward below 1.5/im. As Figure]^ shows, a, z ^ 2.2 galaxy with H/3 -h [O III] lines of combined 
flux of ^ 5 X 10 ^^erg/s/cm^, with a rest-frame Equivalent Width (EW) of ~ 3000A, reproduces 
the observed fluxes, with [O III] 5007A redshifted just red-ward of F140W, contributing the extra 
observed flux to F160W. Similarly, a z ^ lA galaxy with equally strong fla emission could also 
possibly reproduce the observations. 

Our multiple nebular emission lines model strongly constrains the redshifts at which we would 
expect a low redshift interloper to be z 1.4 ± 0.05 (strong Her emission) and 2 ; 2.2 ± 0.05 

(strong H/3 -h [O III] emission). As shown in Figure 2, in order to produce the observe F140W flux 
levels, fainter emission lines are required. This places strong constraints on the wavelength range 
at which strong emission lines are expected. Even if our detailed treatment of nebular emission was 
inaccurate, the observed photometric break implies the existence of some sort of bright emission 
line between 1.2/im to 1.7/im (the combined wavelength coverage of the F140W and F160W filters). 
This is well within the bandpass of the G141 observations. 

As demonstrated here, it is possible to create a scenario in which emission lines of a relatively 
low-z galaxy can mimic the photometric break of a higher redshift object. While these models 
would have to be cherry-picked and might appear to be unphysical and contrived, the simplest and 
most direct way to rule out that any of these scenarios are taking place is by demonstrating the 
absence of strong emission lines in spectroscopic observations of this source. 


4. Slitless Observations 


The confirmation of a low-redshift emission line interloper would imply the existence of a rather 
peculiar galaxy. The HST WFC3 grism mode, and its resolution of a few tens of angstroms per 
resolution element is extremely adept at detecting emission line objects ( Stranghn et al.pOll ) and 
provides us with an opportunity to test the scenario proposed in the previous section. 


As shown in the middle panel of Figure the WFC3 Grism field is extremely crowded and 
slitless spectroscopic observations are therefore difficult. Most spectra, especially the ones from 
fainter sources such as the three lensed images of MACS0647-JD, are contaminated by the spectra 
of nearby objects. On the other hand, the WFC3 G141 grism is well understood, well calibrated, 
and WFC3 has has been remarkably stable over the past six years. We are therefore able to generate 
very accurate two-dimensional simulations of what this field would look when observed with the 
WFC3 G141 grism. The inputs of these simulations are mosaics of the entire MACSJ0647.7-h7015 
field using data from the CLASH program. Mosaics in the F105W, F125W, F140W and F160W 
can then be used to detect the footprint and estimate the spectral energy distribution of every 
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object in the field. Simulations are then generated by using the aXeSIM package, along with the 
photometric data and the latest WFC3 G141 calibration products. With the exception of the still 
poorly calibrated zeroth spectral order, the resulting simulations are a very good match to our 
observations. These simulations are similar to what is done when computing contamination and 
optimal extraction weigths as part of routine grism extraction (e.g. Pirzkal et al.||2004 ; Brammer et 


[ar]|2012 ). The observations and the simulations are shown in the middle and right panels of Figure 
respectively. 

Figure [^qualitatively compares our shallow observations to our simulated dispersed image of 
the strong emission line source already shown in Figure In these simulations, and since these 
three lensed images are partially resolved, we have conservatively set the size of the emission line 
region to be equal to the measured sizes of JDl, JD2 and JDS in the broad band images. These 
sizes were measured from the available broad band imaging and identical to the ones listed in |Coe 


et al. (2013). Our simulated emission lines are therefore more spatially diluted than they would be 


if we had simply assumed an unresolved point source. It is immediately obvious that any emission 
line(s) bright enough to reproduce the observed broad band photometric break would be detected 
in 3 orbit grism observations. As Figure [^ shows, there is no sign of such an emission line in our 
slitless observations for either JDl, JD2 or JD3, demonstrating with little doubt that the CLASH 
z ^ 11 candidate is devoid of any bright emission lines in the wavelength range covered by the G141 
grism. This significantly strengthens the case that this source is a 2 ; 11 Lyman-break galaxy. 


For a more quantitative comparison, we show the extracted one-dimensional spectra of the 
three MACSJ0647-JD sources in Figure As this Figure shows, the low spectral resolution of the 
G141 grism would make the emission lines shown in Figure less prominent and broader. These 
lines would nevertheless remain readily detectable at the > 5a level. As we described in Section 
there are actually narrow ranges of redshifts {Az ^ 0.05) over which our models can reproduce 
the broad-band photometric break. While this would result in the emission line in Figure to be 
shifted, it would remain easily detectable in at least JDl and JD3, where contamination is either at 
longer wavelength, or vertically offset from where we would expect an emission line to be located. 

Figure shows the spectra in units of e~/s with the corresponding flux limits shown on the 
right axis. We show these extracted spectra in these units as their noise level is sky dominated 
and approximately constant 0.02e“/5). Figures shows the extracted spectra of the three 
lensed images as well as that of the LBG presented in Figure [T} The expected flux levels for JDl, 
shown in the top panel of this Figure are 0.02e“/s and the relatively high level of contamination 
(red dashed line) renders the detection of the continuum impossible in these observations. The 
forthcoming observations, chosen to avoid contamination, will triple the integration time and greatly 
improve our odds of detecting the continuum emission from this source. 
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5. Conclusion 


We have presented results based on our shallow, first epoch slitless spectroscopic observations of 
the galaxy cluster MACSJ0647.7+7015. These data have allowed us to obtain the first spectroscopic 
observations of the three lensed images of the z ^ 11 LBG candidate MACSJ0647-JD of|Coe et al. 


(2013). Following a detailed modeling of the entire field of view, we showed that we can exclude 
a scenario in which the photometric break at 1.5/im in the SED of this object could be caused by 
a large EW emission line. Any emission line strong enough to produce the observed broad band 
photometric break would have been easily detected in these observations. We conclude that this 
object remains an extremely strong candidate for a genuine z ^ 11 galaxy. Based on the simulations 
in this paper, we expect the continuum level red-ward of the Lyman Break to be detected in planned 
deeper future observations. 


We would like to thank the referee for his or her comments and suggestions to improve the 
quality of this paper. Support for AZ was provided by NASA through Hubble Fellowship grant 
#HST-HF2-51334.001-A awarded by STScI. This work was supported in part by grant HST-GO- 
13317.13 from the Space Telescope Science Institute, which is operated by AURA under NASA 
contract NAS5-26555. 
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Fig. 1.— The observed broad band flux levels for JDl (error bars) and the spectrum of a 2 ; = 11.2, 
IO^^Mq Lyman Break Galaxy. The model fluxes (blue circles) are in good agreement with the 
observations. 
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Fig. 2.— The observed broad band flux levels for JDl (error bars) and the spectra of two young 
emission line objects at z=2.2 (Top Panel) and z=1.4 (Bottom Panel). The model fluxes (solid 
circles) are in good agreement with the observations, although the observed break is sharper than 
what the models allow when including the nebular continuum emission, as we did here. Models 
that reproduce the observed break can be generated by allowing a strong z ^ lA Ha or a strong 
z ^ 2.2 [OIII] emission line to be the dominant source of the light detected using the HST F140W 
and F160W filters. Several such models at z ^ 1.4 ± 0.05 and z 2.2 ± 0.05 are in principle 
possible. 
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Fig. 3.— Left Panel: Composite (F105W+F125W+F140W+F160W) image of CLASH 
MACSJ0647. The three lensed images of MACSJ0647-JD are shown. Right Panel: The same 
field, dispersed using the HST G141 grism. The locations of the JDl, JD2 and JDS lensed images 
are shown in the left Panel. The locations of the spectra of JDl, JD2 and JDS in the G141 data 
are shown in the right Panel. 
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Fig. 4.— G141 observation of the lensed images JDl, JD2 and JDS. Panels (a) show the model of 
the spectral contamination from all of the sources in the field. Panels (b) show our observations. 
Panels (c) show our G141 simulations of what the lensed images would look like dispersed if its 
spectrum was that of the emission line object shown in FigureAs these simulations oi 3. z ^ 2.2 
H/3+[OIII] emission line galaxy show, an emission line at 1.5/im bright enough to reproduce the 
observed photometric break would be easily detected in our observations of JDl and JD2. We see no 
evidence of such a strong emission line in Panels (d) where we show our contamination subtracted 
observations of JDl, JD2 and JDS. The contamination estimate is accurate and subtracts cleanly 
from the actual observations. The main residual are the two point like zeroth orders at 1.13/im 
and 1.6/im. This is expected as zeroth orders are relatively poorly calibrated and hence only 
approximately simulated. These are however located far away, particularly in the y-direction, from 
where we would potentially expect an emission line from JDl. 
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Fig. 5.— Extracted spectra of JDl, JD2 and JDS using a narrow extraction window of 1.5 
pixel to maximize the signal-to-noise (line with error bars). We also show the emission line model 
spectrum (solid line) of a z 2.2 H/3+[OIII] line interloper in the top panel for JDl. This is the 
same spectrum shown in Figure The expected contamination level is shown using a dash line. 
The contamination is caused by the overlapping spectra of other sources in the field. The emission 
line models for JD2 and JDS are similar to the one shown in the top panel for JDl when accounting 
for the slight difference in spatial size and lower fluxes. 
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Fig. 6.— Extracted spectra of JDl, JD2 and JDS using a narrow extraction window of 1.5 
pixel to maximize the signal-to-noise (line with error bars). We also show the emission line model 
spectrum (solid line) of a high-redshfit LEG in the top panel for JDl. This is the same spectrum 
shown in Figure The expected contamination level is shown using a dash line. 




















